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Two-loop electroweak corrections at high energies 

A. Denner®', M. Melles'^ and S. Pozzorini'^ 

’^Paul Scherrer Institut (PSI), CH-5232 Villigen, Switzerland 

^’Institut fiir Theoretische Teilchenphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 

We discuss two-loop leading and angular-dependent next-to-leading logarithmic electroweak virtual corrections 
to arbitrary processes at energies above the electroweak scale. The relevant Feynman diagrams involving soft- 
collinear gauge bosons 7, Z, W^ have been evaluated in eikonal approximation. We present results obtained from 
the analytic evaluation of massive loop integrals. To isolate mass singularities we used the Sudakov method and 
alternatively the sector decomposition method in the Feynman-parameter representation. 


1. Introduction 

The main task of future colliders such as the 
LHC or an e“''e“ Linear Collider (LC) is the in¬ 
vestigation of the origin of electroweak symmetry 
breaking and the exploration of the limits of the 
Electroweak Standard Model. In order to dis¬ 
entangle effects of physics beyond the Standard 
Model, the inclusion of QCD and electroweak ra¬ 
diative corrections into the theoretical predictions 
is crucial. 

In the energy range of future colliders, i.e. at 
energies above the electroweak scale, ^/s ^ Mw, 
the electroweak corrections are enhanced by large 
logarithmic contributions [I] of the type 

The leading logarithms (LL), also known as Su¬ 
dakov logarithms [2], correspond to M = 2N, the 
next-to-leading logarithms (NLL) to M = 2N—1, 
etc. The logarithmic dependence on various kine¬ 
matic invariants r = s,t,u,... gives rise to sub¬ 
leading logarithms that involve ratios of invari¬ 
ants, and which we denote as angular-dependent 
logarithms 

We will consider the kinematical region |r| ^ 
M^, where all invariants are much larger than 
the electroweak scale. The general form of elec¬ 
troweak logarithmic corrections is complicated by 


the hierarchy of mass scales nit ~ Mh ~ Mz ~ 
Mw ^ mf ^ with heavy masses at the elec¬ 
troweak scale, light-fermion masses m/, and the 
photon mass A as infrared regulator. As a conse¬ 
quence all logarithms of the large ratios Mw /m / 
and Mw/A have to be taken into account. 

All the above logarithmic terms constitute the 
singular part of the corrections in the massless 
limit. They result either as mass singularities 
from soft/collinear emission of virtual or real par¬ 
ticles off initial or final-state particles, or as rem¬ 
nant of ultraviolet singularities after parameter 
renormalization. 

At the one-loop level, it has been proven that 
for processes that are not mass-suppressed at high 
energies the electroweak logarithms are universal, 
and general results have been given [3] and ap¬ 
plied to gauge-boson pair production at the LHC 
[4]. These results are in agreement with vari¬ 
ous explicit diagrammatic calculations for many 
2^2 scattering processes [5-9]. The approxi¬ 
mate size of the one-loop electroweak LL and NLL 
for a typical 2^2 cross section is 
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at ^/s = ITeV, with 1 — 

The LL and NLL have similar size and opposite 
sign resulting in large cancellations. 

Assuming that at high energies the symmetric 
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phase of the electroweak theory can be used, re¬ 
summations of the electroweak logarithms have 
been proposed based on techniques and results 
known from QCD. Fadin et al. [10] have re¬ 
summed the LL by means of the infrared evo¬ 
lution equation. Kiihn et al. have applied results 
from QCD to resum the logarithmic corrections 
to massless 4-fermion processes, e+e“ ^ ff up to 
the NLL [11] and even to the NNLL [12]. It was 
found that at 1 TeV there is no clear hierarchy be¬ 
tween LL, NLL and NNLL, and that the angular- 
dependent logarithms are important. Melles has 
proposed a resummation of the NLL for arbitrary 
processes [13], which relies on the prescription 
of matching a symmetric SU(2) x U(l) theory 
with QED at the electroweak scale. Recently also 
an extension of this resummation to the angular- 
dependent NLL has been proposed [14]. 

All these resummations amount to exponentia¬ 
tions of the electroweak logarithms. The approx¬ 
imate size of the resulting two-loop LL and NLL 
for typical 2 —> 2 processes at -^/s = 1 TeV is 
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and it is clear that in view of the precision ob¬ 
jectives of a LC below the per-cent level these 
two-loop logarithms must be under control. 

All the above resummation prescriptions result 
from matching a symmetric SU(2) x U(l) the¬ 
ory and QED at the electroweak scale, assuming 
that other effects related to spontaneous symme¬ 
try breaking may be neglected at high energies. 
This assumption needs to be checked by explicit 
diagrammatic two-loop calculations based on the 
electroweak Lagrangian, where all nontrivial ef¬ 
fects related to spontaneous symmetry breaking 
are taken into account, in particular (i) the large 
gap Mw ~ Mz A in the gauge sector, (ii) 
the mixing between the gauge-group eigenstates 
B, resulting into the mass eigenstates 7 , Z, 
and (iii) the presence of longitudinal gauge bosons 
as physical asymptotic states. 

The resummation of the two-loop LL has been 
checked for the massless fermionic form factor in 


Refs. [15,16] and for arbitrary processes in the 
massive Coulomb gauge in Ref. [17]. The resum¬ 
mation of the NLL has so far not been confirmed 
by explicit electroweak two-loop calculations. 

A subset of the NLL is furnished by the 
angular-dependent logarithms of type ( 2 ) with 
M = 2N, L = 1. These contributions are 
numerically important [ 8 , 12 ]. At one-loop or¬ 
der, in the t’ Hooft-Feynman gauge, they result 
only from diagrams where a gauge boson is ex¬ 
changed between two external lines. Similarly, 
the angular-dependent NLL at two-loop order can 
be traced back to a relatively small set of Feyn¬ 
man diagrams. This allows us to perform a di¬ 
agrammatic calculation of the two-loop angular- 
dependent NLL for arbitrary processes. The rele¬ 
vant massive two-loop integrals have been eval¬ 
uated in eikonal approximation, and the loga¬ 
rithms have been obtained analytically using two 
independent methods: the first one goes back to 
Sudakov [2], the other one uses sector decom¬ 
position of Feynman-parameter integrals [18-20]. 
A detailed description of this calculation can be 
found in Ref. [21]. Here we summarize the main 
ingredients and results. 

2. Feynman diagrams 

In the following we consider electroweak 
processes^ (p„) ^ 0 , involving n 

arbitrary mass-eigenstate particles. The kinemat- 
ical invariants are denoted as rki = {pk + Pi)'^ — 
2pkPi, and the matrix elements as 

X = (5) 

We restrict ourselves to matrix elements that are 
not mass-suppressed at high energies. In this case 
global gauge invariance implies 

Y^^MIlk)=0i^—j M, a = 7 ,Z,W±,( 6 ) 

where the gauge couplings J“(fc) act as (trans¬ 
posed) matrices on the external-legs k of the ma¬ 
trix element^. 

^As a convention, all particles ipif, and their momenta pk 
are assumed to be incoming. The corresponding 2 —^ n —2 
processes are easily obtained by crossing symmetry. 
^Details concerning our notation can be found in Ref. [3]. 
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Figure 1. Two-loop diagrams with soft-collinear gauge bosons a,b,c = 7 , Z,W* exchanged between 
external legs j^k,l,m = 1 ,..., n. 


In the t’ Hooft-Feynman gauge, the leading 
mass singularities originate from diagrams with 
soft-collinear virtual gauge bosons coupling to ex¬ 
ternal particles. The relevant two-loop diagrams 
are depicted in Fig. 1, where the soft-collinear 
gauge bosons are exchanged between two, three, 
or four of the n on-shell external legs. 

Each loop integral has to be evaluated for all 
different mass assignments that occur in the elec- 
troweak model. For the internal lines we have 
the cases a,b,c = 7 , Z,W^ and for the external 
masses we assume Mw ^ rriext ^ A. 

The Feynman diagrams are evaluated in 
eikonal approximation, i.e. by neglecting mass 
terms and the momenta of the soft gauge bosons 
everywhere in the numerators, apart from the mo¬ 
menta in the couplings of three soft gauge bosons. 
In the massless limit, longitudinal gauge bosons 
have to be substituted by the corresponding 
would-be Goldstone Bosons using the Goldstone- 
Boson equivalence theorem. 


3. Loop integrals in logarithmic approxi¬ 
mation 

In the evaluation of the loop integrals we in¬ 
clude the LL and the angular-dependent NLL, 
and we use mt — Mh — Mz — Mw, i.e. we ne¬ 
glect logarithms of ratios of heavy masses. Two 
analytical methods have been used: the Sudakov 
method [ 2 ] and sector decomposition [18-20], 


which permits to factorize overlapping ultravio¬ 
let or mass singularities in Feynman-parameter 
integrals. Here we only sketch the main steps of 
the second method applied to a generic two-loop 
massive integral. 

Step 1 : The integral is written in Feynman 
parametrization and the denominator is split into 
polynomials according to the hierarchy of scales 
s » r A-P A^ in the diagram^ 


I 

D{x) 


sPsix) + rPr{x) -I- ... -I- PP\{x). 


(7) 


These polynomials have various zeros of the form 

m 

P{x) = ^ XiPi{x) = 0, ( 8 ) 

i=l 

at xi = ... = Xm = 0 , which give rise to mass 
singularities. 

Step 2: In order to separate overlapping sin¬ 
gularities, we decompose the sector [ 0 , 1 ]™ into 
m subsectors Uj with xj > xt^j, and in each 
subsector Vlj we perform variable transformations 
Xi Xjx'^, which remap flj [ 0 , 1 ]"* and permit 

®To extract the angular-dependent logarithms log(s/r) 
with r = we compute the integrals in the eu¬ 

clidean region in various limits of the type s 3> t = n, 
s = t ^ u, etc., where we separate the energy scales in 
various ways. 
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to factorize the variable Xj in 


P{x) 




(9) 


Step 3: Recursive application of step 2 per¬ 
mits to factorize all zeros at all scales, until the 
denominator assumes the form 


D{x) 


sPs (x) Xfc -I- rPr (x) 
k 




+...+A 2 A(x)}n 

m 


where all Feynman parameters that give rise to 
mass singularities are factorized and the polyno¬ 
mials P are non-vanishing. This allows for a sim¬ 
ple power counting of the logarithmically diver¬ 
gent integrations. 

Step 4 : In leading-logarithmic approximation 
the polynomials P can be treated as constants 
P{x) ~ .P(O), and all logarithms of ratios of scales 
can be extracted by analytical integration of the 
singular parameters. 


4. Results 


As a basis for the presentation of our two-loop 
results, we recall the one-loop results for LL and 
angular-dependent NLL given in Ref. [3]. 


4.1. One-loop results 

At one-loop level we have 

A4i = A4o(l + few)- (11) 

The most symmetric form to write the elec- 
troweak logarithmic corrections consists in split¬ 
ting them into few = i^sew + i^sem, with 

few = fewIv^Mw ’ 

fern = few-fewlA=Mw- (^^) 

The symmetric electroweak part (sew) corre¬ 

sponds to the case when the photon mass A equals 
the electroweak scale and reads 


fc=i 






(13) 


+E E fe(fc)fe(oiog—log 


l^k a= 7 ,Z,Wi 




where C®'" = F^/(4c^) -I- represents 

the electroweak Casimir operator. The remaining 
part is a subtracted electromagnetic (sem) contri¬ 
bution originating from the fact that A <C Mw, 


fc=i 




2 log ^ log 

mf. 
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^g(fc)Q(/)log^log-^w 
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4.2. Two-loop results 

Detailed results for the individual two-loop di¬ 
agrams in Fig. 1 are presented in Ref. [21]. These 
diagrams have to be combined as follows 


SM2 = Y^ 

a^b 



c 


+ E 




E pab 

^ jklm 


taking into account all sums over virtual gauge 
bosons a,b,c = 7 , Z, and external legs^ 
j,k,l,m = l,...,n, with appropriate symmetry 
factors. These sums can be simplified by means 
of ( 6 ), and it turns out that the result corresponds 
to the exponentiation of the one-loop corrections 
(13),(14) in the form 


M 2 = Mo exp (few) exp (fem), 


(15) 


where the symmetric electroweak and the sub¬ 
tracted electromagnetic parts exponentiate sepa¬ 
rately. 


5. Discussion 

Our result confirms the exponentiation of the 
electroweak LL obtained with the infrared evolu¬ 
tion equation [10] and in the Coulomb gauge [17]. 
The subset of diagrams (A-C) has been evaluated 
also in the special case of massless external parti¬ 
cles and agreement has been found with the form- 
factor calculation of Refs. [15,16]. The expo¬ 
nentiation of the electroweak angular-dependent 

^In the sums only the contributions from different external 
j 7 ^ have to be considered. 
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NLL agrees with the results of Ref. [12] for mass¬ 
less fermionic processes, and Ref. [14] for arbi¬ 
trary processes. These results were obtained us¬ 
ing matching conditions at the electroweak scale. 

In the following we discuss the idea of matching 
by applying it to our result. If we set A = Mw or 
s = M^, we obtain 

A = <^sem ~ ~ exp [(5gew] *^05 

S = ^ (5sew = 0, M 2 = exp [(5sem] Mo, (16) 


i.e. we observe exponentiation within a symmet¬ 
ric SU(2) X U(l) theory (A = Myy) and QED 
(s = M^). This provides a simple consistency 
check of our result. However, we stress that the 
matching conditions (16) are not sufficient in or¬ 
der to determine the interference terms (Jsewi^sem 
between electroweak and electromagnetic contri¬ 
butions, which vanish at both matching points. 
These two-loop terms are fixed by our full elec¬ 
troweak calculation (s ^ A^) and are 

crucial in order to predict the ordering of the two 
exponentials in (15), which is non-trivial at the 
level of the angular-dependent NLL, since 


[^sew, ^sem] — ^ 




■ ( 17 ) 


The fact that Ssem appears on the right-hand side 
in (15) means that the subtracted electromag¬ 
netic corrections are only sensitive to the elec¬ 
tromagnetic charges of the external particles. 


REFERENCES 

1. M. Kuroda, G. Moultaka and D. Schild- 
knecht, Nucl. Phys. B 350 (1991) 25; 

G. Degrassi and A. Sirlin, Phys. Rev. D 46 
(1992) 3104; 

A. Denner, S. Dittmaier and R. Schuster, 
Nucl. Phys. B 452 (1995) 80; 

A. Denner, S. Dittmaier and T. Hahn, Phys. 
Rev. D 56 (1997) 117; 

A. Denner and T. Hahn, Nucl. Phys. B 525 
(1998) 27. 

2. V. V. Sudakov, Sov. Phys. JETP 3 (1956) 65 
[Zh. Eksp. Teor. Fiz. 30 (1956) 87]. 

3. A. Denner and S. Pozzorini, Eur. Phys. J. C 
18 (2001) 461; Eur. Phys. J. C 21 (2001) 63; 


S. Pozzorini, Dissertation, Universitat Zurich, 
2001, hep-ph/0201077. 

4. E. Accomando, A. Denner and S. Pozzorini, 
Phys. Rev. D 65 (2002) 073003. 

5. W. Beenakker et al., Nucl. Phys. B 410 
(1993) 245; Phys. Lett. B 317 (1993) 622. 

6. P. Ciafaloni and D. Comelli, Phys. Lett. B 
446 (1999) 278. 

7. M. Beccaria et al., Phys. Rev. D 61 (2000) 
073005; Phys. Rev. D 61 (2000) 011301; 
Phys. Rev. D 64 (2001) 053016; Phys. Rev. 
D 65 (2002) 093007; 

M. Beccaria, F. M. Renard and G. Verzeg- 
nassi, Phys. Rev. D 63 (2001) 053013; Phys. 
Rev. D 63 (2001) 095010; hep-ph/0203254. 

8. M. Beccaria, F. M. Renard and G. Verzeg- 
nassi, Phys. Rev. D 64 (2001) 073008. 

9. J. Layssac and F. M. Renard, Phys. Rev. D 
64 (2001) 053018; 

G. J. Gounaris, J. Layssac and F. M. Renard, 
hep-ph/0207273. 

10. V. S. Fadin et al., Phys. Rev. D 61 (2000) 
094002. 

11. J. H. Kiihn, A. A. Penin and V. A. Smirnov, 
Eur. Phys. J. C 17 (2000) 97. 

12. J. H. Kiihn, S. Moch, A. A. Penin and 
V. A. Smirnov, Nucl. Phys. B 616 (2001) 286. 

13. M. Melles, Phys. Rev. D 63 (2001) 034003; 
Phys. Rev. D 64 (2001) 014011; Phys. Rev. 
D 64 (2001) 054003; hep-ph/0104232. 

14. M. Melles, Eur. Phys. J. G 24 (2002) 193. 

15. M. Melles, Phys. Lett. B 495 (2000) 81. 

16. M. Hori, H. Kawamura and J. Kodaira, Phys. 
Lett. B 491 (2000) 275. 

17. W. Beenakker and A. Werthenbach, Phys. 
Lett. B 489 (2000) 148; Nucl. Phys. B 630 
(2002) 3. 

18. K. Hepp, Gommun. Math. Phys. 2 (1966) 301. 

19. M. Roth and A. Denner, Nucl. Phys. B 479 
(1996) 495. 

20. T. Binoth and G. Heinrich, Nucl. Phys. B 585 
(2000) 741. 

21. A. Denner, M. Melles and S. Pozzorini, in 
preparation. 



